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ABSTRACT 



0{N) vector sigma models possessing catastrophes in their action are studied. 
Coupling the limit N oo with an appropriate scaling behaviour of the coupling 
I constants, the partition function develops a singular factor. This is a generalized 

■ Airy function in the case of spacetime dimension zero and the partition function 

of a scalar field theory for positive spacetime dimension. 



1. Introduction 



1^ Matrix models have been understood as representing stochastic triangulated sur- 

O ■ faces and thus interpreted as quantum gravity theories. They are treated in the 

■ "double scaling limit" N ^ oo, g ^ Qc^. The same kind of limits can be applied to 

O(A^) vector models i'iS which are connected with statistical ensembles of polymers. 



These models have been studied with the usual expansion and renormalization 
i! erouD techniaues i'ifl'i. This aoDroach involves a c 



Qh! group techniques H'Iffl'i. This approach involves a certain amount of guess work. 

I If an action possesses a singularity (synonymous: a catastrophe) the double scaling 

limit is naturally defined by letting tend to infinity and the deformation of the 
singularity go to zero. The partition function develops then a singular factor which 
is a generalized Airy function if spacetime has dimension zero or a partition function 
of a scalar field theory if this dimension is nonzero (positive). The Airy function 
depends in general on u scale invariant parameters, where u is the codimension of the 
singularity. 

Singularities can be classified i and form s-dimensional families. The s parameters 
of these families are called "moduli" . If s = the families are discrete and are grouped 
into A, D, and E series. They are related with the corresponding Lie algebras by 
their symmetry. The A series can be realized in one-vector models, D and E series 
need two O(A^) vectors at least. The class of quasihomogeneous singularities is the 
biggest class we have studied so far and shows some hitherto unobserved structures. 

By application of diffeomorphisms singularities can be brought to a canonical form. 
These diffeomorphisms "reparametrize" the coupling constants. The canonical form 
contains the full information of the singularity and defines the universality class of the 
multicritical behaviour. The notation of universality class thus obtains mathematical 



rigour. But the problem of ill-defined (unstable) partition functions can also be 
approached by the theory of singularites. Stability is either obtained by using complex 
contours (as used for the Airy function Bi((^)) or by an appropriate infinitesimal 
deformation. Moreover, renormalization group equations and beta functions can be 
derived a posteriori. 

In a zero dimensional sigma model not all field degrees of freedom participate in 
the singularity. The residual degrees behave Gaussian at the saddle point and are 
integrated over. They give rise to the regular part of the partition function and possess 
a standard expansion. In sigma models with nonzero spacetime dimension small 
momenta are also treated as deformation parameters. Large momenta are included 
in the Gaussian variables. The elimination of the Gaussian variables is a nontrivial 
procedure. It ought to be correct in the singular variables at least to the order of the 
singularity (the method used by Di Vecchia et. al. in IS is wrong). In section 4 we 
describe this elimination in a nontrivial case. 

2. Elementary catastrophes in zerodimensional spacetime 

We outline here the results of 0. Let $ be an O(A^) vector field and be the 
Lebesgue measure on Hat. Then we define the partition function 



Formal, analytic or finite power series are permitted. After integration of the angles 
we obtain 





where g{^) possesses the expansion 
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log X - ^ 



(4) 



We normalize gi to one. 
If we require 



0, k>m+l 



(5) 



we have 




(6) 



and a catastrophe of the type can be produced at Xq if 



/^(xo) = for all 1 < A; < m. 



(7) 



This is achieved if the couphng constants and xq are chosen as the critical values 

Xq = m. (9) 

From (1) we have > only if m is odd. Thus the partition function (|l|) at the 
critical point is ill-defined. However, if instead of @ we have 

gm+i = e > 0, gk = 0, k>m + 2 (10) 

convergence is guaranteed. 

The critical values (^) define the canonical form of the Am catastrophe. We 
define deformations by 

gk = gt + Qk {2<k<m) (11) 
xo = xI + Qq (12) 

so that instead of (^ we have (for 1 < A; < m — 1) 

f^'\xl + Go) = -{k - 1)! E (^^ _ ^jm'-'Qi + 0^(9). (13) 

In fact, translation invariance permits us to submit the deformations (|15), (p^ ) to the 
constraint 

/(™)(x^ + eo) = (14) 

which fixes 0o in terms of {0^}™- 
Now we scale x as 

X — Xq = Xrj (15) 

so that N ^ oo and A — >■ combine to render the (m + l)st order term finite and 
normalized 

2 (m + 1)!'-' ^ ^' m + 1 ^ ^ 

This implies 

2 \ ^^J+T 



If each linear combination (|T^) scales as 



/ 2 \ 



for N —>■ oo then for 1 < < m — 1 



lim ^^/('^Hm + eo) = a (19) 

Af^oo 2 kl 



is finite, provided 

ak = l- — ^ (20) 
m + 1 



and 



h 

771 

Ck = (21) 



The resulting partition function is to leading order in 

N_ 1 

TT 2 N rt \ f 1 \ '" + 1 

Zn = ^Z) e^^^"^^ (]^) • >^(Ci, C2, . . . , U-i) (22) 
where F is the A^-type generalized Airy function 

/■ f m-l ^m+l 

nci,C2,... ,Cm-i)= I dvexp\j2Ckv' + i-ir^—T\- (23) 

■^c(™) [ ^ m + 1 J 

If m is odd, C"-™^ is the real axis with positive orientation. If m is even, IS 
chosen as a linear combination of complex contours each one running from infinity to 
infinity. The singular part of the free energy is 

F,(C)=l0gy(Cl,C2,... ,Cn^-l)• (24) 

It satisfies Airy function differential equations and a renormalization group equation 

where (|T^, ( ]TP| ) has been substituted for C^. The {Ca:} are the "double scale invariant" 
variables. The beta functions {Pk} in (^) can be calculated to the first order in the 



3. Quasihomogeneous singularities 

The singularities of Section 2 involve only one variable. If we want singularities 
with r variables we have to define a corresponding class of more complex models. We 
shall see that models with r vector fields have the desired properties. 

Let 

$2, ...,$.) = i I E 9r{^') + E %((*^*i)') 1 (26) 

[l<i<r l<i<j<r- J 

and 

Zn = J nrf<f.exp{-iV(7($i,$2,... ,$r)}. (27) 

1=1 



We introduce the shorthands 





= $2 
I 






ei 








tij 






Sn- 


_i define the 



V 

Fit) = / n n ^((e. ■ ek? - t,u) 

•' i=l j<k 



(28) 
(29) 

(30) 



(31) 



where only those [j, k) are included in the product for which Sjk is nonzero. If all Sjk 
are nonzero and is big enough, F{t) ( pT]) contains the factor 



61,62, ... ,er 



(32) 



where F*^*"^ is the Gram determinant of the arguments listed. It is a polynomial of 
1 

{tji.} of degree r. For the partition function ( p7[ ) we obtain 



4^ = cpj i[dt,kh{t) 



{n+y 



— exp 



i=i ^* 



(33) 



where h is independent of A^, ID takes account of the Riemann sheets necessary, and 
the phase function (reduced action) ^ is 



log< {W_Xi)T^''\ei,e2, ... , 6^ 



1=1 



l<i<j<r 



(34) 



l<i<r 



The singular saddle point has to be found in ID x (IR^ 
Notation is simplified if we use 



Then the saddle point S has to satisfy 



tjk c'F('-) 



F(^) dt, 



jk 



Xigi\xi) + ^ CijSij{,C,ij) = 1, (for all i) 
CjkSjk'iCjk) = 0, (for all (j, k), j <k). 



(35) 

(36) 
(37) 



In the generic case this saddle point is Gaussian but we are interested in the 
case in which the Hessian has nonzero corank which is achieved by an adjustment of 
couphng constants. With 



— T 



Ui (38) 



dxidxj 
we obtain for i ^ j 

and for i,j,k pairwise different 



/ 2 



a,, (t^j) (39) 



Thus a corank r of the Hessian can be achieved if in addition to (BHl 



Ui = (alH) (42) 
cr.fc = (all(j,fc), j<A;). (43) 

The location of the singularity S is at {xf'\ tjjj} and we use normalized variables 

(0) 

(44) 



For r = 2 we have performed the calculation explicitly I. 

After integration of the Gaussian degrees of freedom {rjk} we have an r variable 
action 

(46) 

and the singular partition function 



2 



(47) 



4i = / d'C h{0 exp 
Assume /(^) is a deformation of a critical function 

/c(0 = /5(0 + />5(0 (48) 



where fs is quasihomogeneous of degree S and type a and fys has higher degree than 
6. To explain these concepts consider a monomial of C, 

e = ei'i2' ■ ■ ■ c'^ (49) 

and apply the multiplicative group 1R+ to ^ by 

6 ^ A^'Cn a^^Q (50) 

so that 

Then is called "of degree 6 and of type a" if 

k-a = S. (52) 

A polynomial is quasihomogeneous if each monomial contained in it satisfies (0). 
Consider as an example the singularity Wio (r = s = 2). Its canonical form is i 

+ (ao + ai6)e?e2' + ^2' (53) 
with ao,ai as moduli (oq 7^ ±2) and the principal quasihomogeneous part 

fs{0=Ct + aoei^l + e2 (54) 

with 

ry=( , 
'4' 6 



a = {],l),S = l (55) 



/>^(0 = «ie?e2- (56) 

We will always adjust a so that 6=1. 
Now we consider a deformation of fdO 

/(0 = EA(0 + />i(0 (57) 

<5=0 

obtained by a change of coupling constants 

{Yn} - {7:^ + ©n = 7n} (58) 

and expansion points 

{6 = O}'; ^ {A,;K. (59) 

We can impose r constraints from translation invariance so that the {Aj}^ are ex- 
pressed linearly in terms of the {6„}. The scaling is done such that the degree 6 = 1 
terms absorb the factor 

— ■ A = 1. (60) 



Each term of lower degree < 5 < 1 



fsiO= E hd' (61) 



k 



can be expanded in the {On} 

V = EA/'a;n0n + O2(e) (62) 

n 

if it is not kept equal to zero by a translational invariance constraint. Then we perform 
the double scaling 

^sr. = (y ) tsr.- (63) 

ve„^o 

In the singular partition function appears the generalized Airy function 

^Z] = [(TieM-Zm (64) 
J c 

m = z,fi. (65) 

0<5<1 k 

(S-k = S) 

Of course the 6 = 1 part of Z{C,) can be normalized to the canonical form. The 
number of arguments of $ is thus reduced to the codimension of the singularity. If 
the sigma model studied is too "narrow", some of these arguments may be missing 
(see Xg singularity in i). The contour C belongs to an equivalence class which runs 
from infinity to infinity and is thus invariant under translations and dilations (pT|). 

In i we derive string equations for these models. The singular partition function 
serves as a ground state functional for the VF-algebra and is invariant under a "tri- 
angular subalgebra" depending only on the type a. If the singularity is homogeneous 
and 

111, 

a = ao{-, , -) (66) 

this subalgebra is the inhomogeneous gl{r) (r translations in coupling constant space 
as invariant subalgebra). 



4. Nonzero spacetime dimensions 

In the case of spacetime dimensions D > we obtain partition functions for r 
scalar fields instead of Airy functions. We will sketch here only the case of the 
singularities (r = 1). A partly incomplete, partly wrong discussion of these cases was 
given by Di Vecchia et. al. in iS. 

Consider the action 

S = j d^x i^-d,^d,^{x) + \p^^\x) + t/(<f ')(x)| (67) 



with the potential 

U{a) = —a". (68) 

n=2 

Then a standard trick gives us the partition function in terms of two scalar fields a 
and p 

Z = J DaDpexp[-NSeff{a,p)] (69) 

with 

Seff = / d^x[U{a){x) — ia{x)p{x)] 



+ ^Tr log [-A + /32 + 2Z/9]. (70) 



The saddle point is assumed at constant values o"o,po 



U'iao) = ipo (71) 



. (27r)^ 

= (3^ + 2ipo. (73) 

The dimensions D are interpolated and firstly restricted to the interval < D < 2, 
so that renormalizations are not necessary. 

The fields a and p fluctuate around their stationary values (0)"(E1) 

a{x) = ao{l + a{x)) (74) 
p{x) = poil + /3ix)). (75) 

In terms of the Fourier transforms of a and /3 the Hessian is 



1 r d^k 
2J (2^ 

-laopo 2plT.{k) j \ p{k) 



(76) 



where 

^« = /(2.)- 
The Hessian is diagonalized by 



m = I S^lip' + m')iip - kf + m')]-\ (77) 



+ X^{k)N^{k)fi{-k)f]{k)} (79) 
where N±{k) are the squares of the eigenvectors 

A^+(A;) = a(A;)^ - 1, N^{k) = 1 - b{kf (80) 

and it turns out that 

a{k)b{k) = 1. (81) 

The eigenvalues X±{k) are 

A±(fc) = i{^"(ao)a2 + 2p^S(fc) 

± [(t/"(ao)ao^-2p^S(A:))2-4aoVg]^}. (82) 

We will make A_(/c) zero at A; = 0. Then X+{k) is positive and ^(fc) a Gaussian 
degree of freedom. If we put the system into a box with periodic boundary conditions, 
then the spectrum of k is discrete and only ?)(0) is non-Gaussian, whereas fi{k), k ^ 0, 
are Gaussian. In this case the saddle point for the Gaussian integration is determined 
by 



di{k) 



0, all k (83) 



g^^O.allMO (84) 

and its location depends on 17(0). Because of translational invariance these equations 
are solved by 

^(A;) = fi{k) =0 for A; ^ (85) 



and with 

by additionally 

which means 



Seff\^^k)=m=0,k^O = Sred{i{0),fl{0)) (86) 

d 



dm 



Sred{m,Vm=0 (87) 



- A+(o)iv+(o)e(o) = y: -^sum,m)- m 

n=3 dm 

Here the superscript n denotes the order in both arguments. 
An iterative solution of (|88D gives 

oo 

m = Him) = T.^ivioy (89) 

1=2 



and 



9. 
n 



oo 



E -^(0)"- (90) 



n=2 



Thus after the Gaussian integration we end up with a zero dimensional treated 



in Section 2. First one calculates the critical coupling constants (7^ in (pOD and in turn 
the critical constants in the potential U{a). 

Now we turn to the case of physical interest, namely infinite volume and a con- 
tinuous momentum spectrum. In the case of the singularity in ( ]90|) the coupling 
constants are deformed as usual 

/n = /^ + e„(2<n<m). (91) 

In addition the domain of small momenta 

|A;|<A (92) 
is considered as deformation, whereas large momenta 

|fc|>A (93) 

belong to Gaussian degrees of freedom ^{k)^f]{k) and are integrated over. Our aim is 
a partition function of a scalar field (j) in conventional form 



D(j)exp 

+ E / d^'x 4>{x)' + ^ / d^x 4>{xr^' I (94) 



which replaces the Airy function (|23|) . Contrary to the Airy function there is no linear 
term in 0, there is a term (^"^, and normalization is applied to the kinetic energy 
and not to 0™+^. Moreover the mass term M^0^ must be positive. This can be 
achieved by approaching the canonical form of the singularity in the neighbourhood 
of deformations on a specific path, namely on a surface of stable Ai singularities 
contingent to A^. 

In the first step we reduce 5'e// to Sredivik)), \k\ < A, that possesses a power 
series expansion analogous to ( PUj) in terms of convolution powers (?})" for n > 3. We 
expand A_(A;) to first order in k"^ and the coupling constant deformations {On} (|9TD- 
Then we scale k by 

k = N-^k' (95) 
X = N^x'. (96) 



The exponent A must be positive to map the domain of small k on ]Ro for large A^. 
As in Section 2 wee firstly examine the term 0™"+^ which entails 



(P{x') = C^"'^N^r]{x). (97) 

With this ansatz we enter the kinetic energy term and find 

m — 1 , , 

2(m+l)-D(m-l) ^ ' 

and 



with 



c(-) = ^2n,(.popi(2-:^)-l (99) 

nn = /^(A:^ + mr". (100) 
The mass is obtained from a double scaling constraint on the coupling constants 

/ \ 1 1 D 

hm E(^-1)0X-' =^fi(2--) — . (101) 

ve—o \n=2 / ^ 

This constraint fixes the limiting contour on the stable Ai surface. 

From (^TJ) we obtain for the intermediate coupling constant Qn in (|90D 

m 

gn-gl = nY^a^:::^Qi. (102) 
1=2 

Then the double scaling limit defines the scale invariant variables Cn 

a = hm (C(-))-"A^>^"'"' f:ai7)e, (103) 

ve—o 1=2 

where the critical indices are 

2(m+l)-D(m-l) ^ ' 

xt^ = 2 A (see (98)). (105) 

The coupling constant Fm+i in ( P^ ) is a critical value and thereby fixed 

F^+i = <7;;+i(C(™))-'"-^ (106) 

However, this coupling constant varies with the mass of the fundamental field $. 



Now we consider the interval of dimensions 2 < D < 4. Then Hi but no n„, n > 2 
( |100| ), diverges. This entails that ctq, po and the bare coupling constants /„ are infinite 
and must be renormalized. First we continue Hi analytically from < D < 2 to 
2 < D < A and call this function 11"" 

Then we replace the saddle point equation (|72D by the renormalized equation 

= Ur = ao-a^ (108) 
^oo = [ i^ipT' (divergent). (109) 



. {27r)D 

We insert ctq (|108|) into the potential (|68D and reorder to obtain 

U'iao) = -Am' + Y.fr«T-' (HO) 

^ n=2 

= ^Am2 + f/'-^"'((T5'^"). (Ill) 

The divergent term Am^ is absorbed in zpg*^" and finally in the renormalized mass. 

After these renormalizations all arguments go through as long as A (pHD is positive, 
i.e. for 

171 -\- 1 

D<D^ = 2 -. (112) 

m — 1 

But this is the classical condition for the partition function (0) to be superrenormal- 
izable. 

For the renormalizable case D = the scaling procedure is not defined due 
to the divergence of the critical indices ( |104|) and A (0). We can however 

renormalize the parameter by introducing 

N\D) = A^2(-+i)-D{™-i) (113) 

so that the scaling procedure defined in terms of N' is analytic in D. Therefore we 
propose to treat the renormalizable case for all 2 < D < 4 by analytic continuation 
in D as well. 

Details about how the results presented in this section are derived and more details 
can be found in i. 
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